V3. Submitted to ApJL 



Estimating Be Star Disk Radii using 
H« Emission Equivalent Widths 

Erika D. Grundstrom and Douglas R. Gies 

Center for High Angular Resolution Astronomy and 
Department of Physics and Astronomy, 
Georgia State University, P. 0. Box 4IO6, Atlanta, GA 30302-4106; 
erika@chara. gsu. edu, gies @chara. gsu. edu 

ABSTRACT 

We present numerical models of the circumstellar disks of Be stars, and we 
describe the resulting synthetic Ha emission lines and maps of the wavelength- 
integrated emission flux projected onto the sky. We demonstrate that there are 
monotonic relationships between the emission line equivalent width and the ratio 
of the angular half- width at half maximum of the projected disk major axis to 
the radius of the star. These relationships depend mainly upon the temperatures 
of the disk and star, the inclination of the disk normal to the line of sight, and 
the adopted outer boundary for the disk radius. We show that the predicted 
Ha disk radii are consistent with those observed directly through long baseline 
interferometry of nearby Be stars (especially once allowance is made for disk 
truncation in binaries and for dilution of the observed Ha equivalent width by 
continuum disk flux in the l^-band). 

Subject headings: stars: emission-line. Be — circumstellar matter — techniques: 
interferometric — techniques: spectroscopic 



1. Ha Spectral and Spatial Emission Models 

Be stars are rapidly rotating B-type stars that lose mass into an equatorial, circumstellar 
disk (Porter & Rivinius 2003). The disk flux is observed in emission lines (especially the 
hydrogen Ha line) , an infrared excess, and the polarization of starlight scattered in the disk. 
It is now possible to investigate the actual extent and geometry of the circumstellar disks 
through optical long baseline interferometry of nearby Be stars (Stee et al. 1995; Quirrenbach 
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et al. 1997; Chesneau et al. 2005; Tycner et al. 2005, 2006). Most of these investigations 
have focused on the angular resolution of the disks in a narrow spectral band centered on the 
Ha emission line. The long baseline interferometric observations are limited to the nearby 
sample of Be stars at present, so we need other indirect methods to help estimate the sizes 
of Be star disks in more distant targets. 

Here we explore how the disk radius may be estimated from measurements of the disk 
Ha emission flux. We begin with an idealized case to illustrate some of the important 
parameters that define the relationship between equivalent width and radius. Suppose that 
the disk is isothermal and geometrically thin and that the line emission is dominated by an 
inner optically thick region that appears as an ellipse projected onto the plane of the sky. 
We also make the simphfying assumption that the specific intensity of the disk emission 7^ is 
isotropic. Then the wavelength integrated Ha flux will depend on the product of the surface 
flux, the wavelength range where Ha is optically thick, and the projected solid angle. 



where < AA > is the wavelength interval over which the line is optically thick (averaged 
over the visible disk), i is the disk normal inclination angle, is the boundary radius of 
the optically thick disk, and D is the distance to the star. We measure this flux through 
the emission equivalent width Wx, which is given in wavelength units relative to the local 
continuum flux. 



where e is the ratio of disk continuum flux to stellar flux in the Ha region, 7| is the stellar 
speciflc intensity near Ha (assumed isotropic), and Rs is the radius of the star. We can then 
equate these two expressions to find the predicted relationship between the ratio of disk to 
stellar radius and Ha emission equivalent width. 



This special case is characterized by a disk radius that varies as the square root of 
the emission equivalent width, and this expectation is supported by the interferometric 
observations of Tycner et al. (2005) who find that the emitted Ha flux scales with the square 

of the disk radius. However, the coefficient of the relationship also depends on the ratio of the 
stellar to disk specific intensities that is dependent upon the stellar and disk temperatures. 
Furthermore, there is an inclination dependence in both the cosi projection factor and the 
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< A A > term (through the range in velocities and Doppler shifts viewed through the disk) . 
The < AA > term also depends upon the assumed disk density since the optical depth is 
based on the product of neutral hydrogen density and the emission line profile. Finally, 
there must also be some emission contribution from the optically thin outer disk that is 
ignored in expression above. Thus, in order to predict accurately the relationship between 
equivalent width and disk radius, we need a more complete model that accounts for the 
specific dependences on disk temperature, inclination, and density. 

Wc have adopted the Be disk model approach of Hummel & Vranckcn (2000) that is 
based upon models developed by Hornc & Marsh (1986) and Hornc (1995) for accretion 
disks in cataclysmic variables. The disk is assumed to be axisymmetric and centered over 
the equator of the underlying star, and the gas density varies as 

p{R, Z) = poR-^'exp 

where it! and Z are the radial and vertical cylindrical coordinates (in units of stellar radii), 
Po is the base density at the stellar equator, n is a radial density exponent, and H{R) is 
the disk vertical scale height. The neutral hydrogen population within the disk is found by 
equating the photoionization and recombination rates (Gies et al. 2006). The disk gas is 
assumed to be isothermal and related to the stellar effective temperature Tes by — 0.6Tefr 
(Carciofi & Bjorkman 2006). 

The numerical model represents the disk by a large grid of azimuthal and radial surface 
elements, and the equation of transfer is solved along a ray through the center of each element 
according to 

/A = 5f(l-e--^) + /^e--^ (5) 

where Ix is the derived specific intensity, iSf is the source function for the disk gas (taken as 
the Planck function for the disk temperature T^), is the specific intensity for a uniform disk 
star (taken as the product of the Planck function for Tcs and an Ha photospheric absorption 
line derived from the grid of Martins et al. 2005), and tx is the integrated optical depth along 
the ray (for assumed Keplerian rotation; see eq. [4] in Home 1995 and eq. [7] in Hummel 
& Vrancken 2000). The first term apphes to all the disk area elements that are unocculted 
by the star (see eq. [12] in Hummel & Vrancken 2000) while the second term apphes to all 
elements that correspond to the projected photospheric disk of the star. The absorption 
line adopted in /| is Doppler shifted according to solid body rotation for the photospheric 
position in a star that is rotating at 90% of the critical value. 

The model code calculates a synthetic line profile over the range from —2000 to +2000 
km s~^ at 10 km s~^ intervals by summing the product of projected area and specific intensity 



2 {h{R)) 



-4- 



over the disk grid. We then integrate the hne profile relative to the continuum to form a 
predicted equivalent width (taken as positive for net absorption or negative for net emission). 
The model also forms a synthetic image of the star plus disk in the plane of the sky by 
summing the intensity over a 2.8 nm band centered on Ha. We collapsed this image along the 
projected major axis to get the summed spatial intensity, and then we determined the radius 
where the summed intensity drops to half its maximum value (excluding the spatial range 
corresponding to photospheric intensity). We adopt this half- maximum summed intensity 
radius as the effective disk radius in the following discussion. 

Wc need to specify the stellar mass and radius in the code in order to determine the 
disk Kcplerian velocities, and we adopted these as a function of stellar effective temperature 
(or spectral subtype) using the eclipsing binary results for B-stars from Harmanec (1988). 
We determined the relationship between model equivalent width and disk radius by running 
a sequence of models defined by the spectral subtype of the star, the disk inclination, the 
disk radial density exponent n, and the selected value for the outer boundary of the disk. 
The final variable for any sequence is the base density po, and we computed models over a 
range in po that corresponds to Ho; profiles with no visible emission up to strong emission 
cases with \Wx\>5Q A. 

We show the predicted relationship between the disk radius and Ha equivalent width in 
Figure 1 for the case of a B2 Ve star (Tcs = 23100 K) assuming n = 3.0 and an outer boundary 
at lOORg. The lower solid line shows the pole-on case (?' = 0°). At higher inclination angles 
{i = 50° and i = 80° for the middle and top lines, respectively) the projected disk area is 
smaller (oc cosi), and consequently the emission equivalent width is smaller for a given disk 
radius. We do not show the i — 90° case because the shear broadening approximation for the 
emission profile width breaks down for an edge-on orientation (Horne & Marsh 1986). The 
symbols along each sequence indicate positions of fiducial base density values. Note that all 
three curves show the expected square root dependence predicted for the simple case (eq. [3] 
above) . 

The three dashed hues in Figure 1 show the corresponding inclination cases for a B8 Ve 
star (Teff = 11600 K) again made assuming n — 3.0 and an outer boundary at 100i?s. These 
curves are all shifted to the left relative to the B2 Ve sequences because the photospheric 
absorption component is larger in cooler B-stars. The model profiles of the B8 Ve sequence 
have less Ha emission at a given disk radius than their B2 Ve counterparts because at this 
lower temperature the disk-to-star intensity ratio, the vertical scale height, and the thermal 
emission line broadening are all smaller. 

The dotted line in Figure 1 shows how the Ha equivalent width is reduced when the 
outer boundary for the disk grid is moved inwards from lOORg to 25Rs for the case of a B2 Ve 
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star with a disk inclination of 80° and a disk density exponent oi n — 3.0. The outer disk 
region contains low brightness, optically thin gas, and although the emission contribution 
from any particular outer radius area element is small, the projected area of such elements 
increases with radius so that the selection of the outer boundary condition is important. 
We found that increasing the outer boundary from lOOi?^ to 200/?^ for this parameter set 
increased the equivalent widths by only 10%, so we adopted the 100i?s boundary as our 
nominal choice. The selection of the boundary radius will mainly be important for Be stars 
in binary systems where the outer disk will be truncated by the gravitational influence of 
the companion. 

Finally, the dot-dashed line in Figure 1 shows another sequence for the same B2 Ve, 
i — 80°, and 100i?s boundary case, but calculated this time with a larger disk density 
exponent of n = 3.5 (i.e., a steeper drop off in density). The overall shape of this curve 
is almost the same as the one for n — 3.0 (although the fiducial base density points have 
very different locations), and this suggests that the particular choice of density law is not 
important for the derived equivalent width - radius relationships. Investigations of the 
infrared flux excess in Be stars suggest that the exponent falls in the range n = 2 — 4 (Porter 
& Rivinius 2003). We caution that the results become more sensitive to the outer disk 
boundary for n < 3 where the emission from the optically thin, outer regions assumes more 
significance. 

We have constructed such model sequences for subtypes BO, B2, B4, B6, and B8 Ve, 
for inclinations i = 0°, 50°, and 80°, and for outer boundary radii of 17Rs, 25Rs, 50Rs, and 
100i?s. These numerical results, an interpolation program, and our code are available to 
interested readers at our Web site^. 



2. Comparison with Disk Radii from Interferometry 

We can test the predictions of the model with the interfcrometric observations of nearby 
Be stars for which the projected disk shape and size in the sky are known. The long baseline 
interfcrometric observations in the narrow band surrounding Ha are fit with an elliptical 
Gaussian, and the reported diameter is the FWHM of the Gaussian fit to the major 
axis. We hst in Table 1 the observed and adopted parameters for the Be stars with Ha 
interfcrometric observations from Quirrenbach et al. (1997) and Tycner et al. (2005, 2006). 
The columns of Table 1 list the star name, a contemporaneous measurement of W\ and a 
code for the data source, adopted effective temperature (Zorec, Fremat, & Cidale 2005), 
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the stellar angular diameter determined from flux fitting methods (Underbill et al. 1979; 
Ochsenbein & Halbwachs 1982), the interferometric disk radius, ratio of the projected minor 
to major axis r, and a code for the data source, an estimate of the disk inclination, and 
finally the ratio Rd/Rs from interferometry and from W\ (plus the adopted outer boundary 
for the calculation of the latter) . 

We determined the ratio Rd/Rs for each of the observations listed in Table 1 by inter- 
polating in the model sequences for T^s, i, and outer disk boundary. We estimated the 
disk inclination angle by equating the predicted value of the projected minor to major axis 
with the observed value. 



where the disk vertical dimension is evaluated at radial distance Rd, Cg is the speed of sound, 
and Vy^ is the Kcplcrian velocity at R — Rg (sec cq. [4] in Hummel & Vrancken 2000). The 
ratio Cs/Vk is ~ 0.022 among main sequence B-stars for the temperatures, masses, and radii 
given by Harmanec (1988). We assumed an outer disk boundary of 100-Rs for all but the 
known binary stars where we adopted the Roche radius of the primary star instead (Gies et 
al. 2006; Hubert et al. 1997). 

The model disk radii predicted from W\ arc listed in Table 1 and compared to the 
interferometric radii in Figure 2. The predicted radii arc generally in reasonable agreement 
with the observed radii, however, the model tends to overestimate the radii for a given 
equivalent width (or to underestimate the equivalent width for a given radius) for the weaker 
emission and cooler stars, rj Tau and /3 CMi. Furthermore, the predicted radii of the stronger 
emission stars appear to be systematically lower than observed. We suspect that this latter 
problem is due to an underestimation of the true emission equivalent width in very active Be 
stars. If Be star disks contribute a fraction of the continuum flux in the V"-band spectrum 
of e = F^/F^, then we need to renormalize the equivalent widths to (1 + e)Wx in order to 
compare them with the model results where the equivalent width is referred to the stellar 
continuum alone. We lack direct information on the value of e for the targets in Table 1, but 
we can approximately estimate the renormalization factor from the work of Dachs, Kiehling, 
& Engels (1988) who found that (1 + e) ^ (1 - O.OOSWx)'^ based upon flts of the spectral 
energy distributions in a sample of Be stars (see their eq. [40]). The revised radii including 
this correction {open squares in Fig. 2) are in better agreement with the interferometric radii. 
Thus, the predictions of the models appear to be generally consistent with the available 
interferometric data. 

The method outlined above offers a new way to estimate disk radius for distant Be 
stars that is based only upon the Ha equivalent width and estimates of spectral subtype, 
inclination, and outer disk boundary. It should prove to be a useful addition to the traditional 
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method from Huang (1972) that is based upon measurements of the separation between 
double peaks of the Ha profile and upon an assumed mass. We caution readers that the 
derived radii for individual targets are probably only accurate to within ±30% at the moment 
(percentage standard deviation of the residuals for the dilution corrected radii in Fig. 2) and 
that higher precision will require detailed modeling of the spectral energy distribution and 
emission lines. Such models will probably need a more detailed physical description of the 
disk (Carciofi & Bjorkman 2006) with attention paid to possible asymmetries in the disk gas 
(Porter & Rivinius 2003). The effort involved is certainly merited in studies of selected Be 
stars, but the simpler approach presented here may be valuable in surveys of Be stars where 
observations are limited to low resolution spectroscopic measurements of the Ha equivalent 
width (or a photometric counterpart; McSwain & Gies 2005). The method will also be useful 
for planning future interferometric observations. 

We are grateful to Dr. Christopher Tycner and an anonymous referee for their insight 
and comments that greatly aided our investigation. This work was supported by the National 
Science Foundation under Grant No. AST-0205297 and AST-0506573. Institutional support 
has been provided from the GSU College of Arts and Sciences and from the Research Program 
Enhancement fund of the Board of Regents of the University System of Georgia, administered 
through the GSU Office of the Vice President for Research. 
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Table 1. Be Stars with Interferometric Ha Measurements 
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Fig. 1. — Plots of the predicted Ha equivalent width and disk half-maximum radius for 
several model sequences. The solid lines show the relationships for models with a B2 Ve 
central star, an outer boundary of lOORg, and inclinations of 0° {bottom), 50° (middle), and 
80° (top). The three dashed hnes show the results from models with the same inclinations 
but with a B8 Ve star. The dotted hne shows how the relation changes for a model with a 
B2 Ve star and disk inchnation of 80° when the outer boundary is reduced to 25Rs. The 
dot-dashed line shows that the relation is almost identical for a model with a B2 Ve star, 
i = 80°, and 100i?s boundary when the disk density exponent is changed from n = 3.0 to 
3.5. Disk base density values of po = 1 x 10^^^, 2 x 10^^^, and 5 x 10^^^ g cm~^ are indicated 
along sequences by a plus sign, diamond, and square, respectively. 
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Fig. 2. — A comparison of the Ha disk radii obtained directly from interferometry with those 
predicted by the models using the Ha equivalent width (diamonds centered on the errors 
bars). The solid line shows the expected one-to-one correspondence. The open squares show 
how the predicted radii increase if the equivalent widths are renormalized to account for 
dilution by the disk continuum F-band flux (Dachs et al. 1988). 



